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Two-nucleon momentum correlation functions are investigated for different single thermal sources 


at given initial temperature (T) and density (p). 


To this end, the evolutions of various single 


excited nuclei at T = 1 — 20 MeV and p = 0.2 - 1.2 po are simulated using the thermal isospin- 
dependent quantum molecular dynamics (T'hIQM D) model. Momentum correlation functions of 
identical proton-pairs (Cpp(q)) or neutron-pairs (Cnn(q)) at small relative momenta are calculated 
by Lednický and Lyuboshitz analytical method. The results illustrate that Cpp(q) and Crn(q) keep 
sensitivities to the source size (A) at lower T or higher p, but almost not at higher T or lower p. And 
the sensitivities become stronger for smaller source. Moreover, the T', p and A dependencies of the 
Gaussian source radii are also extracted by fitting the two-proton momentum correlation functions, 
and the results are consistent with the above conclusions. 


PACS numbers: 
I. INTRODUCTION 


Properties of nuclear matter is one of the most in- 
teresting topics in heavy-ion physics [1-5] and lots of 
works have been done around zero temperature, includ- 
ing the nuclear equation of state (EOS). However, the 
studies on properties of nuclear matter at finite tem- 
peratures are relatively limited. Many previous works 
mainly focus on the temperature dependence of hot nu- 
clear matter and the nuclear liquid-gas phase transition 
(LGPT) [6-15], the ratio between shear viscosity to en- 
tropy density (7/s) [16-20], as well as the nuclear gi- 
ant dipole resonance [21-23] etc. Among above works, 
the relationship between the phase transition tempera- 
ture and the source size has been investigated [6]. In 
Ref. [6], the finite-size scaling effects on nuclear liquid-gas 
phase transition probes are investigated by studying de- 
excitation processes of the thermal sources by the isospin- 
dependent quantum molecular dynamics model (IQMD). 
It explores several probes including the total multiplicity 
derivative, second moment parameter, intermediate mass 
fragment multiplicity, Fisher’s power-law exponent as 
well as Ma's nuclear Zipf’s law exponent, and obtains 
the phase transition temperatures from them. Recently, 
the deep neural network has also been used to determine 
the nuclear liquid gas phase transition [24] and estimate 
the temperature of de-excited nuclei by the charge mul- 
tiplicity distribution [25]. The latter work proposes the 
charge multiplicity distribution can be used as a ther- 
mometer of heavy-ion collisions. 


Considering that the intermediate-state at high tem- 
perature and density in the evolution process of nu- 
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clear reactions cannot be directly measured, one always 
explore properties of nuclear matter and the dynami- 
cal description of heavy-ion collisions through the anal- 
ysis of the final-state products. As well known, the 
two-particle momentum correlation function, or called 
Hanbury-Brown Twiss (H BT) technique in the final- 
state has been used mostly as a probe of the space- 
time properties and characteristics of the emission source 
[26-28]. The two-proton momentum correlation func- 
tion has been explored systematically by a lot of ex- 
periments and different models, several reviews can be 
found in Refs. [29-32]. In various HBT studies, effects 
of the impact parameter, the total momentum of nucleon 
pairs, the isospin of the emission source, the nuclear sym- 
metry energy, the nuclear equation of state (EOS) as 
well as the in-medium nucleon-nucleon cross section has 
been discussed in the literature [33-39]. Even more, nu- 
clear structure effects were also carefully investigated, 
such as the effects from binding energy and separation 
energy of the nucleus [40], density distribution of va- 
lence neutrons in neutron-rich nuclei [41], as well as high 
momentum tail of the nucleon-momentum distribution 
[42] etc. Two-proton momentum correlation function 
was also constructed in few-body reactions as well as a- 
clustered nucleus induced collisions [43-47]. In addition, 
two light charged particles momentum correlation func- 
tion also offers a unique tool to investigate dynamical 
expansion of the reaction zone [39]. 


Here we extend the HBT method of final-state in- 
teraction to study the time-spatial information of the 
finite-temperature nuclear system. The purpose of the 
present paper is to systematically investigate the re- 
lationship between two-particle momentum correlation 
functions and system parameters, such as the source- 
temperature, density as well as system-size in a frame- 
work of the thermal isospin-dependent quantum molec- 


ular dynamics (ThIQMD) model [6, 15, 18]. In addi- 
tion, Gaussian source radii are quantitatively extracted 
by assumption of Gaussian source fits to the HBT dis- 
tributions. In this article, the evolution process of the 
excited nuclear source at an initial temperature varying 
from 1 MeV to 20 MeV are studied. The present work 
selects six different systems with similar ratio of neutron 
to proton numbers, i.e., N/Z ~ 1.3, which include (A, Z) 
— (36,15), (52,24), (80,33), (100,45), (112,50), and 
(129,54). Then, Lednickj-Lyuboshitz theoretical ap- 
proach [48] is applied for calculating two-particle momen- 
tum correlation functions constructed based on phase- 
space information from the evolution process of the single 
excited nuclear source by the T'AIQM D model. 

The rest of this article is organized as follows. In Sec- 
tion IT, we firstly describe the thermal isospin-dependent 
quantum molecular dynamics model (TAIQM D) [15, 
18], then briefly introduce the Hanbury-Brown and Twiss 
(H BT) technique using Lednickÿ and Lyuboshitz ana- 
lytical formalism. In Section 111, we show the results 
of the ThIQMD plus the LL method for the source- 
temperature dependence of two-particle momentum cor- 
relation function. The two-particle momentum correla- 
tion function and the influences of different temperature 
and system-size are systematically discussed. The de- 
tailed analysis of the extracted Gauss source radii are 
presented under different source-temperature and den- 
sity. Furthermore, the momentum correlation function of 
two-neutron is also analyzed. Finally, Section IV gives a 
summary of the paper. 


IIl. MODELS AND FORMALISM 
A. THE ThIQMD MODEL 


Since the two-particle momentum correlation functions 
are calculated in a specific application of the Lednicky 
and Lyuboshitz theoretical simulation, the true single- 
particle phase-space distribution at the freeze-out stage 
is required. In this paper, the thermal isospin-dependent 
Quantum Molecular Dynamics transport model is used 
as the event generator, which has been applied success- 
fully to study the LGPT [6, 25]. In the following discus- 
sion, we introduce the modified model briefly. As well 
known, isospin-dependent Quantum Molecular Dynam- 
ics (IQM D) model was used to describe the collision 
process between two nuclei. The Quantum Molecular 
Dynamics transport model is a n-body transport theory, 
which describes heavy-ion reaction dynamics from inter- 
mediate to relativistic energies [49-53]. In the present 
work, we use a single excited source in the ThIQMD 
which is different from the traditional /QM D. Usually, 
the ground state of the initial nucleus is considered to 
be T = 0 MeV in the traditional /QM D model. How- 
ever, the T'AIQM D model described by Fang, Ma, Zhou 
and Liu in Ref. [6, 15, 18], is used to simulate the single 
nuclear source at different temperature T > 0 MeV. 


The main parts of QM D transport model include the 
following issues: the initialization of the projectile and 
the target, nucleon propagation under the effective po- 
tential, the collisions between the nucleons in the nuclear 
medium and the Pauli blocking effect. In the T'AIQM D, 
instead of using the Fermi-Dirac distribution for T' = 0 
MeV, Pi(r) = A [3r?p:(r)] 1/3 the initial momentum of 
the nucleons is sampled by the Fermi-Dirac distribution 
at finite temperature: 
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where the kinetic energy ex = 4—, p and m is the mo- 
mentum and mass of the nucleon, respectively. g (ex) = 
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In addition, the chemical potential u is determined by 


the following equation: 
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In the ThIQM D model, the interaction potential is 
also represented by the form as follows: 


U = Usky + Ucou + Uyug + Usym +Umpr, (3) 


where Usxky, UCoul, UYuk, U Sym» and Ujyp,; are the 
density-dependent Skyrme potential, the Coulomb po- 
tential, the surface Yukawa potential, the isospin asym- 
metry potential, and the momentum-dependent interac- 
tion, respectively. Among these potential, the Skyrme 
potential, the Coulomb potential and the momentum- 
dependent interaction can be written as follows: 

Usi, = AE) + 8, (4) 
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where p and p, are total nucleon density and its normal 
value at the ground state, i.e., 0.16 fm *, respectively. 
The above parameters a, B, and y with an incompress- 
ibility parameter K are related to the nuclear equation 
of state [54-60]. 
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where p, and py are neutron and proton densities, re- 
spectively, 7, is the z th component of the isospin degree 
of freedom for the nucleon, which equals 1 or —1 for a 
neutron or proton, respectively. And Csym is the sym- 
metry energy coefficient. Uc; is the Coulomb potential 
where V. is the Coulomb potential for protons. 


Umpr = ô: n? (e- (Ap)? +1): m (7) 
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where Ap is the relative momentum, 6 and € can be found 
in Refs. [49, 50]. Their value of the above potential pa- 
rameters are all listed in Table I: 


TABLE I: The value of the interaction potential parameters. 


a B y K ô € 
(MeV) (MeV) (MeV) (MeV) ((GeV/c) ?) 
—390.1 320.3 1.14 200 1.57 500 


B. LEDNICKY AND LYUBOSHITZ 
ANALYTICAL FORMALISM 


Next, we briefly review the method for the two-particle 
momentum correlation function proposed by Lednicky 
and Lyuboshitz [48, 61, 62]. The main framework of 
the HBT technique is based on the principle as fol- 
lows: when they are emitted at small relative momenta, 
the two-particle momentum correlation is determined 
by the space-time characteristics of the production pro- 
cesses owing to the effects of quantum statistics (QS) and 
final-state interactions (FSI) [63]. Therefore, the two- 
particle momentum correlation function can be expressed 
through a square of the symmetrizied Bethe-Salpeter am- 
plitude averaging over the four coordinates of the emis- 
sion particles and the total spin of the two-particle sys- 
tem, which represents the continuous spectrum of the 
two-particle state. 

In this theoretical approach, the final-state interactions 
of the particle pairs is assumed independent in the pro- 
duction process. According to the conditions in Ref. [64], 
the correlation function of two particles can be written 
as the expression: 
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where r* = x1 —x» is the relative distance of the two par- 
ticles at their kinetic freeze-out, k" is half of the relative 
momentum between two particles, S (r*, k*) is the proba- 
bility to emit a particle pair with given r* and k*, i.e., the 
source emission function, and Vy (r*) is Bethe-Salpeter 
amplitude which can be approximated by the outer so- 
lution of the scattering problem [65]. In the above limit, 
the asymptotic solution of the wave function of the two 
charged particles approximately takes the expression: 


Wye (r*) = e^ / A (x 
G (p, à) 
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In the above equation, 6, = arg I (1 + iA) is the Coulomb 
s-wave phase shift with À = (k*ae) ^ where a, is the two- 
particle Bohr radius, A, (A) = 2zA [exp (21A) — 1] * is 
the Coulomb penetration factor, and its positive (neg- 
ative) value corresponds to the repulsion (attraction). 
G (p, A) = y Ae (A) [Go (p, A) + iFo (p, À)] is a combina- 
tion of regular (Fo) and singular (Go) s-wave Coulomb 
functions [61, 62]. F(—iA,1,i£) = 1 + (—iA) (i£) /AP + 


(—iA) (—idA + 1) (i£)? /2!2 +... is the confluent hyperge- 
ometric function with € = k*r* + p, p = k*r*. 


—1 


fe (kY) = | Ke (*) — 2h (A) — dk* A, QU) (10) 
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is the s-wave scattering amplitude renormalizied by 
the long-range Coulomb interaction, with h(A) = 
X2 EZ [n (n? + d2)]~" —C—In [A] where C = 0.5772 is 
the Euler constant. Ke (k*) = & + idok” +Pk* +--> is 
the effective range function, where do is the effective ra- 
dius of the strong interaction, fo is the scattering length 
and P is the shape parameter. The parameters of the 
effective range function are important parameters char- 
acterizing the essential properties of the FSI, and can 
be extracted from the correlation function measured ex- 
perimentally [39, 65-67]. 

For n—n correlation functions which include uncharged 
particle, only the short-range particle interaction works. 
For p— p correlation functions, both the Coulomb interac- 
tion and the short-range particle interaction dominated 
by the s-wave interaction are considered. 


Ill. ANALYSIS AND DISCUSSION 


Within the framework of the thermal isospin- 
dependent quantum molecular dynamics model {6, 
15, 18], the two-particle momentum correlation func- 
tions are calculated by using the phase-space infor- 
mation from the freeze-out stage of the excited nu- 
clear source at an initial temperature varying from 1 
MeV to 20 MeV and/or density varying from p = 
0.2po to 1.299. This work performs calculations for 
thermal source systems with different mass including 
(A,Z) = (36,15), (52, 24), (80, 33), (100, 45), (112, 50), 
and (129, 54). 

We firstly calculated the proton-proton momentum 
correlation function Cpp(q) for a finite-size systems at 
temperatures ranging from 1 to 20 MeV. In Fig. 1, the 
results of C,,(q) for temperature of 2, 4, 6, 8, 10 and 
12 MeV at different values of density (0.2po - 1.209) 
are presented. The proton-proton momentum correla- 
tion function exhibits a peak at relative momentum q — 
20 MeV/c, which is due to the strong final-state s-wave 
attraction together with the suppression at lower relative 
momentum as a result of Coulomb repulsion and the anti- 
symmetrization wave function between two protons. The 
shape of the two-proton momentum correlation functions 
is consistent with previous experimental data in heavy- 
ion collisions [68]. For protons which are emitted from 
the lower temperature (T « 8 MeV) source in Fig. 1 (a)- 
(c), the general trend is very similar. The figure shows 
that Cpp(q) increases as p increases for fixed T (T < 8 
MeV). The increase of the density indicates that the 
geometrical size becomes smaller for a source with fixed 
neutrons and protons, which makes the strength of the 
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FIG. 1: The proton-proton momentum correlation functions 
(Cpp(q)) at different densities (i.e., 0.200, 0.4po, 0.60, 0.8po, 
1.0po, and 1.299) for the smaller nucleus (A—36, Z—15) with 
fixed source-temperatures T = 2 MeV (a), 4 MeV (b), 6 
MeV (c), 8 MeV (d), 10 MeV (e) and 12 MeV (f), respec- 
tively. The freeze-out time is taken to be 200 f m/c. 
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FIG. 2: Similar to Fig. 1, but at different source- 
temperatures (T = 2, 4, 6, 8, 10 and 12 MeV) with different 
fixed densities, namely p = 0.2p0 (a), 0.4 po (b), 0.6 po (c), 
0.8 po (d), 1.0 po (e), and 1.2po (f). 


momentum correlation function larger. Finally, the p- 
p momentum correlation function become almost one at 
q > 60 MeV/c. For larger T (T > 8 MeV) in Fig. 1 
(d)-(f), the difference of Cpp(q) between different densi- 
ties becomes smaller. From Fig. 1, it is found that the 
Cpp(q) is almost the same after T = 8 MeV for different 
densities and the p — p momentum correlation function 
becomes almost one above approximately q — 30 MeV/c. 
It indicates that the emitted proton is not affected by the 
change of density when the source temperature beyond 
certain value (T ~ 8 MeV in present work). In order 
to understand which one of the two factors (1.e., temper- 
ature and density) has larger influence, the two-particle 
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FIG. 3: The neutron-neutron (n-n) momentum correlation 
function (Cnn(q)) in the same conditions as Fig. 1. 
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FIG. 4: Same to Fig. 1, but for a larger system 
(A—129,Z—54). 


momentum correlation in fig. 2 is plotted by exchang- 
ing of the two input parameters. From fig. 2, we can 
intuitively observe dependence of the two-particle mo- 
mentum correlation on the source temperature. The de- 
pendence of C, (q) on the source temperature is stronger 
than on density. In other words, the Cpp(q) is more sensi- 
tive to T' than to density p. In addition, for larger p from 
fig. 2 (a) to (f), the difference of Cy, (q) between differ- 
ent densities becomes bigger. Next, we explore whether 
the phenomenon exists in momentum correlation of the 
uncharged-particle pairs. Fig. 3 presents the neutron- 
neutron momentum correlation functions (Chn(q)) for 
temperature of 2, 4, 6, 8, 10 and 12 MeV at different 
values of density, respectively. For neutron-neutron mo- 
mentum correlation function, it peaks at q ~ 0 MeV/c 
caused by the s-wave attraction. Although the Cnn(q) 
has different shape compared with the p — p momentum 
correlation, it has the similar dependence on the source 
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FIG. 5: 


Cpp(q) of different source size system at fixed temperatures (i.e., from left column to right column, they correspond 


to T = 2, 4 and 6 MeV, respectively) or fixed densities (i.e., from top row to bottom row, they correspond to p = 0.2po, 0.4po, 


0.6po, respectively). 


temperature and density. The similar trend in C,,(q) 
and Cyn(q) shows the close emission mechanism in the 
evolution process. 


Fig. 4 shows the results of a larger system at differ- 
ent source-temperature and density and a similar behav- 
ior of C5, (q) is demonstrated. We also observe that the 
proton-proton momentum correlation in larger-size sys- 
tem ((A, Z) — (129,54)) in Fig. 4 becomes weaker com- 
pared to the smaller-size source ((A, Z) = (36,15)) in 
Fig. 1. In view of the above phenomenon, Fig. 5 de- 
scribes the relationship between system-size and momen- 
tum correlation function in more details. The decreasing 
of C, (q) as the system-size increasing for a fixed value of 
T or p can be clearly seen in Fig. 5 (g), which is consistent 
with the previous results of Gaussian source [38, 39, 69]. 
In Fig. 5 (a)-(i), with larger temperature or lower den- 
sity, the difference of C,,(q) between different T or p 
becomes smaller, respectively. The Gaussian source radii 
are extracted for further discussion later in this article. 


From the above plots, we can extract Cinaz(q), i.e., 
the maximum value of C,,(q) as well as the full width at 


half maximum (FWHM) of Cy, (q) distribution, i.e. at 
Cpp(q) = [Cmas(q) — 1]/2. The source-temperature T de- 
pendence of Cmaz(q) and FWHM for the proton-proton 
momentum correlation function with different density are 
given in Fig. 6. As shown in Fig. 6 (a) and (b), both 
Cmaz(g) and FWHM decrease gradually with increas- 
ing T. In addition, both of them increase gradually 
with increasing density. At high temperature, the change 
of Cinaz(q) and FWHM is very small and not plotted 
in the figure. Of course, the behavior of the C. (q) 
and FWHM with T and p can also be cleanly seen in 
Fig. 2, and the increasing of Cmaz(q) and FWHM are 
generally inversely proportional to Gaussian radius ro as 
shown later. Similarly, the system-size À dependence 
of Cmaz(q) and FWHM for the proton-proton momen- 
tum correlation function at T = 2MeV and p = 0.6po 
is shown in Fig. 7. The dependence of C,,,,(q) and 
FWHM on system-size A is quite similar to the temper- 
ature dependence in Fig. 6. The Cmaz(q) and FWHM 
values become smaller for systems with larger size. 


Fig. 8 shows the source-temperature, density, and 
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FIG. 6: Source-temperature T dependence of Cmaz(q) (a) 
and of FWHM (b) of C5, (q) distributions at different densi- 
ties (0.2po - 1.2p0) for the (A = 35, Z = 16) system. 
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FIG. 7: Cmaz(q) (a) and FW HM (b) for different source- 
size systems at given T' — 2 MeV and p — 0.6po. 
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FIG. 8: Gaussian source radius as a function of temperature 
at different density (p = 0.2po, 0.490, 0.6po, 0.8po, 1.0po, 
1.290) for a fixed source size. Panels (a) and (b) correspond 
to the smaller source size with (A — 36, Z — 15) and the 
larger source size with (A — 129, Z — 54), respectively. 


system-size dependence of Gaussian radii extracted from 
two-particle momentum correlation functions, where 
panels (a) and (b) are results with the smaller source 
size and the larger source size, respectively. The radii 
are extracted by a Gaussian source assumption, i.e., 
S(r) ~ exp[-7?/(4r2)], where ro is the Gaussian source 
radius from the proton-proton momentum correlation 
functions. The theoretical calculations for Cpp(q) was 
performed by using the Lednickÿ and Lyuboshitz an- 
alytical method. The best fitting radius is judged by 
finding the minimum of the reduced chi-square between 
the TAIQM D calculations and the Gaussian source as- 
sumption. It can be observed that the Gaussian radius 
increases initially and then decreases slightly after taking 
a maximum value at a certain temperature in Fig. 8. As 
the density decreases, the decreasing speed of the Gaus- 
sian radius of the small system is larger than that of 
the larger system. Fig. 9 shows the Gaussian radius of 
the different system-size varies with the temperature in 
panels (a)-(c) or density in panels (d)-(f). The Gaussian 
source radius is consistent with the system-size, i.e., at 
higher temperature or larger density, the differences of 
Gaussian source between different system sizes are big- 
ger in the low density and low temperature region, but 
the difference in opposite conditions almost disappear. 


From the above discussion, it is demonstrated that 
the strength of the two-particle momentum correlation 
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FIG. 9: Gaussian source radius as a function of temperature 
or density at different source-size systems. Left and right 
column correspond to the density (p = 0.290, 0.6p0, 1.0p0) 
and the temperature (T = 2, 6, 10 MeV), respectively. 


function is affected by the source temperature, density 
and system size. The two-particle momentum correlation 
function strength is larger for a single source with lower 
temperature, higher density or smaller mass number as 
shown in Fig. 1-5. Otherwise, the strength becomes 
smaller. To some extents, the strong correlation between 
two particles is mainly caused by the closing position of 
each other in phase space in both coordinate and momen- 
tum. Varying only one in the three condition parameters 
(temperature, density and system size), lower temper- 
ature means smaller momentum space, higher density 
means smaller coordinate space and small system size 
also mean smaller coordinate space to keep fixed den- 
sity compared with large system size. The dependen- 
cies of the two-particle momentum correlation function 
strength on the source temperature, density and system 
size could be explained by the change of the phase space 
sizes. Two particles emitted from small phase space will 
have strong correlation and those from large phase space 
will have weak correlation. For example, the increase of 
the C,,(q) strength with the increase of the density for 
a fixed system size could be explained by the decreasing 
of the coordinate space as shown in Fig. 1 (a). And the 
small C, (q) strength at temperature higher than 8 MeV 
could be caused by the large momentum space compared 


with lower temperatures as shown in Fig. 1 (d-f). The 
decrease of the Cpp(q) strength with the increase of the 
system size for a fixed density could also be explained 
by the increasing of the coordinate space as shown in 
Fig. 5 (g). Thus it could be concluded that the phase 
space size for the emitted nucleons have strong effect on 
strength of the two-particle momentum correlation func- 
tion, which can also be seen in the extracted Gaussian 
radii as shown in Fig. 8. However, the sensitivities to 
the source radii seem to be different in the different re- 
gions of temperatures and densities. For example, the 
temperature sensitivity are better in the lower region at 
fixed density, and so does for density sensitivity at lower 
temperatures. 


IV. SUMMARY 


In summary, the two-particle momentum correlation 
functions for a single excited source are investigated using 
the Lednicky and Lyuboshitz analytical formalism with 
the phase-space points at the freeze-out stage for differ- 
ent initial temperature and density in a framework of the 
ThIQMD transport approach. We mainly performed a 
series of studies focused on the varied effects of source 
temperature, density and system-size on the two-particle 
momentum correlation functions. The results reflect that 
the shape of the two-proton momentum correlation func- 
tion is in accordance with the previous experimental data 
in heavy-ion collisions [68]. At the same time, the trend 
of the relationship between the two-proton momentum 
correlation and system-size is consistent with previous 
simulation [38, 39, 69]. At low source-temperature, the 
larger density makes the two-particle momentum corre- 
lation stronger. However, at higher source temperature, 
the effect becomes almost disappear. Both protons and 
neutrons have the similar behavior. This work also shows 
that the emitted particles are not influenced by density 
above a given temperature in a single excited source. In 
the same way, the emitted particles are softly influenced 
by temperature below a given density in single excited 
source. In one word, the dependence of the two-particle 
momentum correlation function on the source tempera- 
ture, density and system size could be explained by the 
change of the coordinate and/or momentum phase space 
sizes. In the end, the Gaussian radius is extracted to 
explore the emission source size in single excited system. 
Gaussian radius becomes larger in the larger system. The 
dependence of the extracted Gaussian radius on source- 
temperature and density is consistent with behavior of 
the two-proton momentum correlation function as dis- 
cussed in the texts. 
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